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ÇSJ i After a short abstract introduction on the time évolution driven by non self-adjoint hamil- 
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tonians, we show how the recently introduced concept of pseudo-fermion can be used in 
the description of damping in finite dimensional quantum Systems, and we compare the 
results deduced adopting the Schrôdinger and the Heisenberg représentations. 



I Introduction 



In a séries of papers, [I]-[H], we have considered two operators a and b, with b ^ a\ acting 
on a Hilbert space "H, and satisfying the commutation rule [a,b] = 1. A nice functional 
structure has been deduced under suitable assumptions, and some connections with physics, 
and in particular with what is usually called quasi-hermitian quantum mechanics and with the 
technique of intertwining operators, have been estabhshed. FoUowing Trifonov, j9], we have 
called pseudo-bosons (PB) the particle-like excitations associated to this structure. A similar 
analysis has also been carried out for what we have called nonlinear pseudo-bosons (NLPB) in 
[ID]-[l2], and most of the original results have been recovered also in this more gênerai situation, 
where the main ingrédient is not the commutation rule between a and b but their raising and 
lowering properties, when applied to two fixed biorthogonal bases. The analytical treatment of 
both PB and NLPB turns out to be particularly difïicult in the case where regularity is lost, 
that is, see below, when thèse bases are not Riesz bases. In this case, in fact, the intertwining 
operators appearing in the game (whose square roots are metric operators in the sensé of the 
literature on quasi-hermitian quantum mechanics, see |T3l [Hl [15] and références therein) turns 
out to be unbounded. For this reason, a large amount of mathematical care is required, and 
this makes the rigorous treatment of the physical System rather complicated, [T2] . 

More recently, [16], we have introduced a similar extension for the canonical anticommu- 
tation relation (CAR), following an original idea by Trifonov and his coUaborators, [T7], also 
briefly considered in |T8] and in [19]. The idea is very similar to the previous one: we consider 
again two operators a and b, with 6 7^ a^, acting on a Hilbert space Ti, and satisfying the anti- 
commutation rule {a, 6} = 1, as well as a"^ = b"^ = 0. Of course, iib = we go back to ordinary 
CAR. The functional structure that we deduce out of thèse operators is nice and friendly, since 
T-L is finite-dimensional. The related excitations are called pseudo-fermions, (PF). One of the 
most important features of PF is that ail the problems arising when dealing with PB, because 
of their (in gênerai) unbounded nature, are now absent. 

In this paper we focus on a particular aspect of PF, i.e. on their possible use in the analysis 
of damped quantum Systems. More in détails: in the description of damping in quantum optics 
certain finite dimensional, non self-adjoint, matrices are quite often used. Thèse matrices are 
usually called effective hamiltonians. Solving the related Schrodinger équation, one sees that 
the wave-function of the system goes to zéro for large t, at least if the parameters of the model 
are chosen properly. Already in [16] we have seen that one such a model, originally proposed 
in [20], can be easily written in terms of PF. However, we have not explored the damping 
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features of that system, because that was not our major interest. Here, on the other hand, we 
focus our attention exactly on this aspect, and we work both in the Schrodinger and in the 
Heisenberg représentations, showing that analogous conclusions can be deduced. We also extend 
our analysis to a higher-dimensional Hilbert space, which we will relate to a two-dimensional 
family of PF. For both thèse examples we deduce damping §. 

The paper is organized as follows: in the next section we briefly discuss some basic facts on 
PF and on the time évolution in the Heisenberg représentation, when this is driven by some 
non self-adjoint hamiltonian. This, we believe, is useful since not ail the readers are familiar 
with this kind of time-evolution, and surely less readers are familiar with pseudo-fermions. 
Section III is devoted to our first, two-dimensional, example, while in Section IV we consider a 
four-dimensional, example. For both thèse examples we compare the well known Schrodinger 
point of view with the modified Heisenberg évolution, and we show that, not surprisingly, the 
same conclusions can be deduced. Particularly relevant for us is the fact that, after a somehow 
standard transformation of the effective hamiltonian, we will naturally driven to introduce 
pseudo-fermions in the game. Section V contains an abstract generalization of the approach, 
while our conclusions are given in Section VI. 



II Pseudo-fermions and dynamics 

We begin this section by briefly reviewing the main définitions and results concerning PF in 
one dimension. The extension to higher dimensions will be discussed later on. The starting 
point is a modification of the CAR {c, c^} = cc^ + c = 1, {c, c} = {c\ c^} = 0, between two 
operators, c and c\ acting on a two-dimensional Hilbert space "H. The CAR are replaced here 
by the foUowing rules: 

{a,b} = l, {a,a} = 0, {b,b} = 0, (2.1) 

where the interesting situation is when b ^ . Thèse rules automatically imply that a non zéro 
vector, Lpo, exists in T-L such that a^po = 0, and that a second non zéro vector, "^q, also exists in 
n such that 6t^o = 0, [E]. 

Under thèse conditions it is possible to recover similar results as those for PB. In particular, 
we first introduce the following non zéro vectors 

ipi := bipo, = a^^i/o, (2.2) 
^It might be worth stressiiig that we are not claiming here that PF are always relevant to explain damping. 
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as well as the non self-adjoint operators 

N = ba, m = N^ = a)b^. (2.3) 



We further introduce the self-adjoint operators S^p and S^^ via their action on a generic f eH: 

1 1 

= Y^^Vn. f) Vn, S^f = Y^i^n, f) ^n- (2.4) 

n=0 n=0 

Hence we get the following results, whose proofs are straightforward and will not be given here: 
1. 

acpi = ipo, b^^^i = ^'o. (2.5) 

2. 

NLpn = mpn, m^Hn = ^i^^n, (2.6) 

for n = 0, 1. 

3. If the normalizations of (fo and \l/o are chosen in such a way that {ipo, \E'o) = 1, then 

{iPk,^n) = h,n, (2.7) 

for k,n = 0, 1. 

4. Sp and S^j, are bounded, strictly positive, self-adjoint, and invertible. They satisfy 

\\SJ<\\ipof+yif, ll^^ll < ll^of + ll^ill', (2.8) 

Sp'^n = ^n, S^(pn = "^n, (2.9) 

for n = 0, 1, as well as S*^ = S^^. Moreover, the following intertwining relations 

S^N = ^S^, Spm = NS^, (2.10) 

hold. 

The above formulas show that (i) and behave as fermionic number operators, having 
eigenvalues and 1; (ii) their related eigenvectors are respectively the vectors of J-",^ = {ipo,(pi} 
and J-"^ = {\l/o,\I/i}; (iii) a and are lowering operators for J^^ and J^is, respectively; (iv) 
h and are rising operators for J^^p and J-"^ respectively; (v) the two sets J-"^ and J-"^ are 
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biortho normal; (vi) the well-behavecQ operators and Sq, maps J-"^ in J^q, and viceversa; (vii) 
and S'^ intertwine between operators which are not self-adjoint, in the very same way as 
they do for PB. 

We refer to [16] for further remarks and conséquences of thèse définitions. In particular, for 
instance, it is shown that J-"^ and are automatically Riesz bases for "H, and the relations 
between fermions and PF are discussed. 

II. 1 The algebraic dynamics 

In [20] a non self-adjoint effective hamiltonian is considered in the description of the time 
évolution of the wave-function of a two-dimensional quantum system. The resulting dynamics 
shows a decay of this wave-function. This is just a single référence where such an approach is 
considered. Other références are, for instance, [T7|, [21] and, for more abstract considérations, 
[22] . The starting point is always the same: a Schrodinger équation 




where H^ff ^ -^e//- Actually, in [20] and [17], i^e/f is assumed to be pseudo-hermitian. We 
will not consider this extra assumption in the first part of this section. It is worth stressing here 
that formula (12.111) should be considered as an assumption in itself, since in standard quantum 
mechanics the Schrodinger équation is usually (if not always) assumed to hold for a self-adjoint 
hamiltonian. Of course, if -f^e// does not dépend explicitly on time, the formai solution of (I2.1ip 
is \l/(t) = e"*^'=-''-''*\l/(0). It is clear that, if -ffe// is a finite-dimensional matrix, this exponential 
can be defined using, for instance, a norm convergent séries: 



On the other hand, if i^e// is an unbounded operator, defining e~*^'=-'-^* is quite a délicate 
task, since we cannot even use the spectral theorem. We will not consider this aspect of the 
theory here, since we are concerned only with finite matrices. Solving (I2.1ip means that we 
are adopting the Schrodinger représentation. To move to the Heisenberg représentation, we 
assume, as it is always done in ordinary quantum mechanics, that the time évolution of the 
mean values of the observables do not dépend on the représentation chosen. In other words, 

^i.e. self-adjoint, bounded, invertible, with bounded inverse 




z^{t) = H^ff^it) 



(2.11) 




fc=0 
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if X is an observable of our physical System, calling Xeff{t) its time évolution driven by Hefj, 
we require that 

{^{t),x^{t)) = {^{o),x,fjm{o)) . 

Since this equality must hold for ail possible choices of ^1/(0), we deduce that 

X,ff{t) := é^^-^ii 'X e^'^^ff *. (2.12) 

This is the time évolution of X in the Heisenberg représentation. Notice that, if Hf,ff = H\^j =: 
H, we go back to the usual formula and to the unitary évolution, X{t) := é^^X e~*^*. Notice 
also that, while {XY){t) = X{t)Y{t), {XY)effit) ^ Xeffit)Yeff{t): the time évolution in fl2:T2D 
is not an automorphism of the set of the observables of the System. Another interesting feature 
of (12.121) is that this time évolution is stable under the adjoint: {Xf,ff{t)y = {X'^)f.ff{t). 

Also, it is clear that the Heisenberg équation of motion for Xeff{t) should be modified. For 
that, let us define the foUowing effective commutator hetween two generic operators A and B: 

[A,B]eff:= AB - A. (2.13) 

The effective commutator satisfies, for example, the following equality: given three generic 
operators A, B and C, then 

[AB, Clff = A[B, C]e// + [A, Clfj B + A(Ct - C)B, 

which reduces to the standard formula when C = C\ but not in gênerai. Also, not surprisingly, 
the Jacobi identity ceases to be true for [., .]eff- Still we have [A,B]ljj- = —[A^,B]eff. 
The differential équation for Xeff{t) is the following: 

Xeffit) = [X^jfit), H^fjl^^ = e</*[X, H^j^lj^e-^^^ff^ (2.14) 

which returns the standard équation if H(,ff = H\j^. This formula shows, among other things, 
that an intégral of motion, Z, is not an operator commuting with -ffe//- What Z has to satisfy 
is clearly the following requirement: [Z, Heff]eff = 0. In fact, if this is true, then Zeffit) = 0. 

In [21] the authors have considered a slightly différent point of view, splitting H^ff into 
two parts, both self-adjoint: Hejf = \ {H^ff + hIj^ ~ ^Yi {j^eff ~ ^eff^ ■ ^^i^ ^^J^ rather 
than introducing an effective commutator, they deduce a differential équation for Xeff{t) in 
which both a commutator and an anti-commutator appear. 
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Something more can be said if Hefj is crypto-hermitian, see [H], i.e. if a positive self- 
adjoint, time-independent, operator exists, bounded with bounded inverse, such that H^ff = 
Q~^Hlf^Q. In this case we get the foUowing: 

and the differential équation is 

I {Q^'X^ffit)) =ie^^^ff'[H^ff,Q-'X]e-^^^ff\ 

Let us remark that here, first of ail, the commutator is the standard one. Also, at a first sight 
the time évolution for Q~^X might appear as an automorphism. Nevertheless, it is not hard 
to check that this is not true, exactly because of the présence of 6~^. In fact, 9^^Xe//(t) 
is, in gênerai, différent from {Q~^X)eff{t). A simple case when they coïncide is when H^jf 
is self-adjoint. The essential reason why formulas above look interesting is because only i^e/f 
appears, and not H\^^. In particular, the differential équation shows that Xf,ff{t) is constant 
when Heff commutes not with X but with Q~^X. This is the formula which replaces, in this 
particular case, the one deduced above, [X, Heff]eff = 0. 

Because of the properties of O, we know that 0^^/^ can be defined, and h := ©^/^//g//©"^^^ 
is a self-adjoint operator: h = . Let us now introduce the following automorphism of Bil-L), 
the set of the bounded operators on l-L: je{x) := 0~^/^a;G^/^. This map is invertible, and it 
is clear that h = ië^(-^e//)- further define the following standard Heisenberg évolution, 

a*{X) := e^^^Xe~^^^, X G BiTi), which has ail the usual properties of the time évolution for 
Systems driven by self-adjoint hamiltonians, we deduce that 

Xeffit) = Je' {a\jeiX))) , (2.15) 

which has a first interesting conséquence: this effective time évolution coïncides with a* if this 
latter commutes with je. Moreover, the interprétation is the following: to find Xeff{t), we have 
first to map X into j@{X). Then we let this operator evolve using h, and finally we transform 
back this resuit, using j'q^; therefore, the effective time évolution appears to be the composition 
of thèse three maps. 
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III An example from the literature 



The first example is taken by [T71 [20] . and its relation with PF was already partly considered 
in [16]. The starting point is the Schrodinger équation 

t4f{t) = H,ff<ilit), with H,ff = U~'y ), (3.1) 

where 7a, 7ft > and f G C. 

III. 1 Schrodinger représentation 

With a simple change of variable = e^*\l'(t), F = |(7a + 76), we deduce that i^{t) = H^{t), 
where 

-(7;| (::;:! 

Here I2 is the two-by-two identity matrix, and 7 = |(7a — 7b). The differential équations for 
$o(0 and $i(t) are easily deduced. Let us introduce fi := jt'p — 7^. Then 



Mt) = -nMt), 
$i(t) = -fi$i(f). 



(3.2) 



If fi = then the functions $o(0 a-nd '^i('^) linear in t, so that 



$i(t) i \ e-(^"+^*H(Ai + Sit) 



where Aq,Ai,Bq and Si are fixed by the initial conditions. In particular, if \l/(0) = 

we find that Aq = ipo, Ai = ipi, Bq = —7 — ivipi and Bi = •y cpi — ivipo. 
The norm of \l/(t) is simply the norm in 7/ = C^. Therefore 

||^(t)f = e-*(7.+7.) (1^^ + Bot\^ + \Ai + , 

and this squared norm goes to zéro when t — )■ 00. 

Let us consider now the case in which fi > 0. In this case the solution can be written as 



^(t) = e 



-ila+'ïb) 
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Ao cos( Vfi t) + Bo sin( Vfi t) 
Al cos{VUt) + Bi sm{VUt) 



where, again, Aq,Ai,Bq and Bi are fixed by the initial conditions: Aq = ipQ, Ai = ipi, Bq = 
{—^ (fQ — i V Lfi) and Bi = {-y ipi — iv lpq). In this case is the product of a 

decaying exponential and an oscillating function. More explicitly, 

||^(t)f = e-*(>+^'') (^IAoCos(v^t) + Sosin(v^t)P + |Aicos(v^t) + 5isin(v^t)|2) , 

which decays to zéro independently of the particular values of lpq and (fi. 

Let us finally consider the case in which fi < 0. In this case the solution can be written as 



vi/(t) 



-(7<i+7b) 



Aoexp{^y\^l\ t) + Boexp{-^y\n\ t) 
Al exp(^/|l7[t) + Bi exp(- Y^|fi[t) 

Here, Aq,Ai^Bq and Bi are related to the initial conditions as follows: 



-HPQ-ivipi)\ , Bq = -\ifQ + ^={'yipQ + iv(p 





^1 . , , ^ - ^ ^ V^o) j , Bi = -\^i + (-7 içi + iv v9o) 

In this case it is not completely clear that 11^1/(^)11 goes to zéro for increasing t. However, to get 
such a conclusion, it is enough to consider the asymptotic behavior of e~^^°-'^'^''^* g2y^t^ which 
is the only dangerous contribution appearing in ||\Ef(t)||^. The conclusion is that 11^1/(^)11 — > 0, 
when t — > oo, if 2y^|TÏ[ < 7a + 7è, which is always satisfied. 

The outcome of this analysis is, therefore, the foUowing: the wave-function \E'(t) solving the 
Schrôdinger équation ( IJ. ij) goes to zéro independently of the value ofQ. 



III. 2 Heisenberg représentation 

It is easy to see that fi is related to the eigenvalues of H, which are X± := ±-^|uP — 7^. Hence 
A+ = A_ = if fi = 0, while X± := i^fi when fi > 0. Finally, X± := ±i \/\fî\ when fi < 0. In 
ail thèse cases the eigenstates of H can be written as 




In fact, Hr]± = X±'r]±. Moreover, in gênerai r]+ and are not orthogonal: (?7+,?7_) = 
^7 — z-y/fij, which is zéro only if 7 = {H = W) or if 7 = i\/Vt [H = —W). Due to 
the relation between H and H^ff we also find that 

H.ffT]^ = E±r]±, E± = -'-{-fa + 7fe) ± Vn. (3.3) 



9 




In [16] we have analyzed, in connection with PF, the following non self-adjoint operator, 

u iô 

which coincides with our H simply taking 6 = 2^ and u = 2v. As in [16], it is therefore possible 
to introduce two operators a and b, such that {a, b} = 1, = 6^ = 0, and H = Q (6a — 1 1). 
Introducing also the number operators N = ba and ^ = N\ and assuming, just to fix the 
ideas, that f2 > 0, we can write H = Q (^N — 1 1) and W = Q (OT — 1 1) . Now, we want to 
recover here the same damping we have found working in Schrodinger représentation. For that, 
it is natural to consider the time évolution of the number operators and which should be 
computed as in (12.121) : 

Due to the relation between H^ff and H, H^ff = H — iFl, and to the above expressions for H 
and i/^, we find that, for each observable X, 

X,ff{t) = e-2r*e*f^^* x e"*^^*. (3.4) 

Therefore, recalling that the anticommutation rules imply that A^^ = A^, = 9T, and so on, 
we find that 

while ^f.ff{i) is simply its adjoint. Then, if we estimate the norm of Nf,ff(t), it is trivial to 
deduce that ||A^e//(t)|| < Se"^^*, which goes to zéro when t diverges. Hence, as expected, we 
recover damping also in Heisenberg picture, with no need of introducing extra assumptions on 
the parameters of the system. More connections with PF are discussed in [T6] . 



IV A two-dimensional example 



We are going now to repeat the same analysis considering the following four-by-four effective 
hamiltonian: 

/ -a/3{ui -0J2) \ 

^ _ i U1-U2 -{a + P){ui - UJ2) 

a — P —U2 au2 (3u2 — aui 

V -l3uj2 l3'^uj2 auj2 - fiuji ) 

where a ^ (5 and uji > UJ2 > G- It is évident that iï^e/f 7^ ^Iff show that, also for this 
hamiltonian, a decay can be deduced. 
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IV. 1 Schrôdinger représentation 

As before, we want to solve first the Schrôdinger équation = Hf.ff'^{t) and check whether 
damping is deduced, and in which conditions. Even if iï^e// has no direct physical interpréta- 
tion, still it is interesting, since, with a simple transformation, we can transform the original 
Schrôdinger équation into a new one, with a traceless hamiltonian H, which is also not self- 
adjoint. This is very close to what was done, for instance, in [T71 [THl EU], and in the previous 
section. For that it is enough to proceed as in the previous example, defining 
H = iTl^ + Heff, and T = 2{a-i3) ('^i ~ '^2)5 where I4 is the identity matrix in "H = 
we get 



^. Then 



where 



H 



a 



(3 



( \{a + P){u:i-u:2) 

-UJ2 



-al3{uji - UJ2) 
i(a + /3)(wi -W2) 
auj2 
















i(a-/3)(a;i+a;2) / 
(4.1) 



and 



$2(t) 

As in our first example, the matrix H is traceless. The related Schrôdinger équation of motion 
produces, first of ail, $o(^) = ^'î'o(Oî where fi := 



U2\. Then <î>o(t) = e^* + e 



with and i?o fixed by the initial conditions. The second component of ^\{t\ can be 
deduced from $o(^) since 



^2) 



- (a + - UJ2)%{Ï) - (a - /3)<î'o(t) 



Moreover 




A2 e-('^i+^2)*/2 + X2 + B2 e 



-m 



(4.2) 



As for the asymptotic behavior of '^{t), we recall that 0Ji> uj2> 0. Of course, to conclude that 



$(t) —7- when t — )■ 00, it is enough to consider the worse contribution, i.e. the one 



11 



coming from $3(t): in fact, if this goes to zéro, then it is clear that ail the other components 
go to zéro as well, so that ||\E'(t)|| — when t — > oo. We have, for t ^ 1, 

\E'3(t) ~ exp I ( ^(wi + ^2) - r ] 1 1 = exp | — ^— {au2 - l3u)i) 



,2 / J [a- P 

Recalling that a > (3, \l/3(t) — ?■ if and only if au2 — (3ui < 0, i.e. if 

In other words, when condition (14. 3p is satisfied, not only "^sit), but itself decreases to 
zéro for t diverging. 



IV. 2 Heisenberg représentation 

Let us now consider the foUowing operators 

p\ 

1 



ai := - 



02 



/ 


/33 


a-/3 


0-/3 


-13 




a-p 


0-/3 


o?-fi 


/3(-a2+/3) 






\ 





/3 




a-/3 




/3 


-P^ 


a(a-/3) 


a{a-lS) 


—a 


0? 


a-P 


0-/3 






a{a-P) 


q(o-/3) 



1 




J 



a—p a— fi 

1 ^ i 

/3) a— /3 a 



\ (9) a— /3 a 





— 





1 \ 



/ ^ ^ ^ \ 

a(a—p) a—p a 

1 ~i i 



- / 

a / 



\ 



a(a—l3) a— {3 a 

1 ^ i 

a(a—l3) a— fi a 

^^00 

a— fi a— fi 



They satisfy the foUowing anticommutation rules: {aj,bk} = (5j,fcl4, with {aj,ak} = {bj,bk} = 
0, j,k = 1, 2§|. Therefore they are two-dimensional pseudo-fermions. H can be written in terms 
of thèse operators as 

u ■ ( h , u (^1 + ^2 ^ 
H = z Wi Oi ai + a;2 02 «2 I4 



As in the previous example, the time évolution of a given operator X, can be written in a way 
which extends (13. 4p : 



Œ^fft 



^(2r+i^l+aJ2)tg(wiîni+W20^2)t ^(LUiNi+uJ2N2)t 



^ Other examples of 4 x 4 matrices satisfying thèse rules can be found in [23] 
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where Nj = bjaj and = Nj , j = 1,2. Recalling that [A^i, A^2] = and that = Nj, we get 
Ni,eff{t) = e-(2r+-2)t + ^^(g-it _ ^ ^ 9T2)(e"^* - 1) + 9^2 A^2(e"^* - 1)) , 

and a similar expression can be deduced for A^2,e//(^)- Moreover, Dîj_e//(t) is simply the adjoint 
of Nj^effit). Using ||iVj|| < 1 and < 1, j = 1,2, we get 

||A^i,e//(t)|| < exp{-(2r -ui- U2)t}, 

which goes to zéro, when t — )■ oo, if and only if 2r > ui + uj2, i.e., if condition (14.30 is 
satisfied: we recover exactly (and not surprisingly) the same conclusion as in the Schrodinger 
représentation. 

Remark:— It might be interesting to notice that the matrices Uj and bj above can be related 
to the foUowing, standard, two-dimensional fermion annihilation operators, 

/0100\ /0010\ 



Al 




1 

V / 



A, 



-1 


V / 



as in [16]: aj = TAjT ^, bj = TA^jT ^ , j = 1,2, where T is the foUowing non-singular matrix: 



T 



/ a /3 \ 
110 
a 1 
\ (3 a / 

Here a ^ and a ^ j3. Therefore, this example can be considered as an explicit realization of 
the équivalence theorem stated in [16], concerning the relations between fermions and pseudo- 
fermions. 



V An abstract generalization 



We will now describe how the examples considered in the previous sections can be further 
extended. 

Let us consider two sets of pseudo-fermionic operators, aj and bj, with bj o)-, j = 
1,2, . . . , N , satisfying 



{aj, bk] = Sj^k 1, {%, ûfc} = {bj, bk} = 0, 
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for ail j and k, and let (/5o,o,...,o and \I^o,o,...,05 with (</?o,o,...,05 ^0,0,. ..,o) = 1, be the vacua of the a/s 
and bj's respectively: 

aj^o,o,...fi = ^î^o,o,...,o = 0, 
for ail j. Then, extending the procédure of Section II, the sets J-"^ and J-"^ of the functions 

,n ,n — /i"l?i"2 A"iV,o 

'^n ■— Vni,n2,...,nN " "2 "A^ V^O,0,...,Oî 

and 

rij = 0, 1, are biorthonormal bases of the 2^-dimensional Hilbert space Un- 

Let now = Hf,ff'^{t), with i^e// 7^ -^e//' original Schrodinger équation. As 

in Sections III and IV, we introduce a new wave function, = e^*\l/(t), where F is a real 
constant, to be fixed. The Schrodinger équation for is 

i^it) = H^^t), Hm := H,ff + zri^, 

where Iat is the identity operator on I-Ln- The examples discussed before show that, in some 
interesting situations, the hamiltonian has the foUowing gênerai form: 

N N 

j=i j=i 

where Nj = bj aj is the j — th pseudo-fermionic number operator, while the Qj^s are, in gênerai, 
complex quantities. Let us now define the foUowing quantity: 

k,n k 

where the last equality is a conséquence of the orthonormality of J-',^ and J-"^, and of the fact 
that (fk is eigenstate of Hjsf. Tjsf appears a natural extension of the trace of an operator to 
the présent settings, where a single orthonormal basis is replaced by two biorthonormal bases. 
From (15. ip it is now easy to check that Tn = for ail = 1, 2, . . .. 

In the (extended) Heisenberg représentation, see Section fil. Il the time évolution Xeff{t) = 
a given operator X, assume the following form: 

N N 

j=i 1=1 
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Focusing now to the time behavior of the number operator Nk, and recalling that Nj = Nj and 
Dî| = DTj, for ail j, we deduce that 

TV N 

{Nk)eff{t) = e-t(2r+3(no+-+3(o^)) -Q (^^^ + (^^ulf- _ ^ j j "Q ^ ^^-^tn, _ ^^^^ 

i=i 1=1 

(5.3) 

Damping can be found by estimating the norm of the varions {Nk)eff{t). For that, we first 
recall that \\Nj\\ = < 1, for ail j. The simplest sitnation is when the complex freqnencies 
are ail real: in fact, if Q'(^^j) = for ail j, we easily conclude that 

ll(iV;^)e//(^)ll<3^^e-2^^ 

which goes to zéro when t diverges for ail possible values of positive F. Hence, in this case, 
damping is deduced. It is probably more interesting to consider the situation in which not ail 
the Qj^s are strictly real. In particular, we consider here the opposite situation, i.e. the one 
in which '^{^j) 7^ for ail j, since this is the most dangerous case. Now, the estimate above 
should be replaced by the foUowing: 

N 

\\{Nk)effm < e-(2r+3(f70+-+3(n^))t -Q (2 + e*^(^^))' , 

j=i 

which, clearly, extends the previous one. In this case, damping is recovered if 

1 ^ 

F>-5^|î3(fi,)l, (5.4) 

i=i 

while a similar (and simpler) condition is recovered if only some of the Qj^s are not purely real. 
Notice that the inequality (15. 4p becomes quite simple, F > 0, when 3(f2j) = 0, for ail j. In this 
case, we recover our previous conclusion. 

If, in analogy with Tn, we introduce the generalizediia.ce also for i^e//, recalling that Tn = 
we find that 

k 

Therefore, in order for H^jf to produce damping, the foUowing inequality should be satisfied 
by the effective hamiltonian i^e//: 

TV 

j=i 
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We end this section with some remarks on the pseudo-fermionic structure related to this 
abstract system. The starting point is the pair of operators S^p and S^^ which extends those of 
(12.41) . They satisfy a multi-dimensional version of (I2.10p . 

Nj = mj , mj Sq, = Nj , 

for ail j. Moreover, since the square root of both and Sq, can be defined, we can also 
introduce the self-adjoint operators rij = S]^"^ Nj S^J"^ = nj. Therefore, if the Qj are real, we 
deduce that S*^ = Ht, S*^ and that 



is a self-adjoint operator. 

It is well known that the operator S]^"^ , as well as its inverse 5"^^, can be used to define a 
différent scalar product in Hat. Then the conclusion is the following: a damping effect in "H^r, 
endowed with its natural scalar product < ., . >, can be obtained using certain non self-adjoint 
effective hamiltonians which, when considered again in "Htv, but endowed with a différent scalar 
product, are self-adjoint. In fact, let us define a new scalar product (., .)^ on l-L^ as foUows: 



for ail f,g G "Hat. Then, since the intertwining relation above implies that Hj^ = Hj^ Sqj, 
we get 



for ail f,gE Un- Hence Hj^ turns out to be self-adjoint with respect to (., .)^. 

Alternatively (and more straightforwardly), formula (15. 5 p shows that the effective hamilta- 
nians considered here and producing damping are similar (but not unitarily équivalent) to self 
adjoint operators. A natural and quite interesting question to consider is therefore: is this a 
gênerai requirement? We have no answer now, but this is part of our works in progress. 

VI Conclusions 

After some gênerai remarks on the time évolution driven by a non self-adjoint hamiltonian, we 
have shown how this can be used in the analysis of the decay of finite dimensional Systems 




(5.5) 





/^S'J'HMf,S'J'g) = {HMf,g) 



16 



in connection with pseudo-fermions. In particular, in the examples discussed here, we have 
seen that the effective hamiltonians describing damping are related, in a very direct way, to 
pseudo-hermitian operators. We have also shown that the Heisenberg représentation can be 
conveniently adopted in the analysis of the time évolutions of thèse Systems. 

The next natural step of our research will be to check how much of this approach can be 
extended to infinité dimensional spaces and to try to answer to the following, rather gênerai, 
question: can any effective damping be described in terms of some pseudo-hermitian hamilto- 
nian? 
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